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Abstract 

We investigate the interplay of quantum interference effects and electronic- vibrational coupling in 
electron transport through single-molecule junctions, employing a nonequilibrium Green's function 
approach. Our findings show that inelastic processes lead, in general, to a quenching of quantum 
interference effects. This quenching is more pronounced for increasing bias voltages and levels of 
vibrational excitation. As a result of this vibrationally induced decoherence, vibrational signatures 
in the transport characteristics of a molecular contact may strongly deviate from a simple Franck- 
Condon picture. This includes signatures in both the resonant and the non-resonant transport 
regime. Moreover, it is shown that local cooling by electron-hole pair creation processes can influ- 
ence the transport characteristics profoundly, giving rise to a significant temperature-dependence 
of the electrical current. 

PACS numbers: 73.23.-b,85.65.+h,71.38.-k 



I. INTRODUCTION 



In the past decades, electronic components have become continuously smaller and reached 
the nanoscale. At these scales, electron transport can no longer be understood on purely 
classical grounds but quantum mechanical effects, such as, for example, the quantization of 
energy levels, quantum interference effects and tunneling processes need to be taken into 
account. Single- molecule junctions, where a molecule is contacted by two macroscopic (in 
most cases, metallic) electrodes, represent prime examples of nanoelectronic devices. These 
junctions allow to investigate molecules under controllable nonequilibrium conditions and 
may facilitate an ultimate step in the miniaturization of nanoelectronic devices [l-4]. It is 
thus highly desirable to understand the electron transport properties of these junctions. 

Experimental studies employing a variety of techniques, such as mechanically controlled 
break junctions |5-14], scanning tunneling microscopy 15-25] and electromigration 
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291 ]. have shown that molecular junctions exhi 



phenomena, including, for example, switching 



arge variety of interesting transport 
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}it a 

negative differential resistance 



39_|, diode- and/or transistor-like behavior 
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4l|. Many of these transport 



properties are not fully understood yet. On one hand, this is due to ambiguities in the contact 
geometry, which may lead to large variations in the recorded current-voltage characteristics 
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42]. While strategies to control the contact geometry exist, for example, by the use 



of special linker groups 



33 



45 



43] , direct binding schemes 22|, |44j , or click-chemistry approaches 



46] . the electrodes themself may show irregular behavior. It is therefore expedient 



to identify transport properties of molecules that are robust with respect to the contact 



geometry of the contact, such as, for example, inelastic electron tunneling spectra 
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50j). On the other hand, transport through the molecular conductor represents a complex 



coupling to numerous 
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52] . In contrast 



many-body problem, where, besides multiple electronic states, the 
vibrational degrees 

to other nanoelectronic devices (e.g. quantum dots 53]), electronic- vibrational coupling is 
typically rather strong in molecular systems. This is due to their small size and mass and 
can lead to pronounced vibrational effects in the respective transport characteristics jllL 
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56] . Thereby, the interplay of electronic and the vibrational 



degrees of freedom of a molecular conductor often involves strong nonequilibrium 



as well as quantum mechanical effects Q, [52|. 
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A fundamental quantum mechanical effect is quantum interference. The occurrence and 
observability of quantum interference effects in electron transport through molecular junc- 
tions has received great attention recently [fj], [ 
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60|. Quantum interference effects in 



molecular junctions arise because, very similar to a double slit experiment, an electron may 
take different pathways through a molecule. Thereby, pathways can be spatially different 
paths in molecules with an appropriate topology and/or originate from different paths in en- 



ergy space, e.g., for molecules that exhibit quasidegenerate electronic states 
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52|, 
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Recently, strong experimental evidence for quantum interference effects in molecular junc- 



tions has been reported 
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591 ] . Moreover, quantum interference effects have been observed 



in the closely related field of electron transport through quantum dots that are set up as 



Aharonov-Bohm interferometers 



6(3-681. 

Depending on the magnetic flux that is threading 



the quantum dot, the transmission amplitudes of the electrons interfere constructively or 



destructive 



y with each other, leading to sizeable oscil 



device 67H69| . A great deal of theoretical work 61 



ations in the conductance of such a 

en devoted to 
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study quantum interference effects in molecular junctions, which is, on one hand, due to 
their fundamental importance but, on the other hand, also due to possible device applica- 
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84J, thermoelectric devices 85|, |86j or spin 



tions, such as, for example, in transistors 

filters 87|, |88|. Thereby, in most of the studies, the effect of electron-electron interactions 
and electronic-vibrational coupling has been neglected. 

In this article, we investigate the influence of electronic-vibrational coupling on quantum 
interference effects in single-molecule junctions. Thereby, we extend our earlier studies 
[ylH] and analyze in detail the mechanism of vibrationally induced decoherence. Moreover, 
we consider the effect of local cooling by electron-hole pair creation processes and show how 
these processes facilitate a mechanism that allows to control quantum interference effects 
by an external parameter, that is, the temperature of the electrodes. Only recently, this 
mechanism has been experimentally verified by Ballmann et al. jl^] for a number of different 
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junctions, including different molecular species. Note that, in contrast to Refs. 
we do not consider the vibrational degrees of freedom as being part of a thermal bath but, 
on the contrary, as active degrees of freedom. 

For our studies, we employ a nonequilibrium Green's function approach developed by 



Galperin et al 
brational 



9jJ , which was recently extended by some of us to account for multiple vi- 



92| as well as multiple electronic degrees of freedom 
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94j. The approach is 
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based on a separation of electronic and vibrational time scales and, therefore, allows a non- 
perturbative description of electronic-vibrational coupling as it is required in this nonequilib- 
rium transport problem. Moreover, it suits to describe quasidegenerate levels, which exhibit 
pronounced quantum interference effects. Besides this approach, a variety of other methods 
has been developed to describe vibrationally coupled electron transport through nanoelec- 
tronic devices. This includes other nonequilibrium Green's function approa ches , wh ich are 



based on eith er pertur bation theory 95Ml02j ] or nonpertur 



Dative schemes 



tering theory 106Ml09| . maste r equatio n methodologies 
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loa - lioa i. scat 



numerically exact schemes 
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116( | and a number of 



The article is organized as follows: In Sec. Ill we introduce the model Hamiltonian 
that we use to describe vibrationally coupled electron transport through a single-molecule 
junction. The nonequilibrium Green's function approach, which we employ to calculate the 
corresponding transport characteristics, is outlined in Sec. Ill Bl In Sec. IIII A\ we introduce 
and analyze in detail basic quantum interference effects that occur in single-molecule junc- 
tions. This facilitates the discussion of the effect of vibrationally induced decoherence in the 
following sections, Sees. IIII Bl — IIII El Thereby, we study first simplified electronic-vibrational 
coupling scenarios that allow to investigate basic decoherence phenomena due to vibrations 
in the resonant (Sec. IIII Bp and the non- resonant transport regime (IIII CI) . More realistic 
models are studied in Sees. IIII Dl and IIII E\ where we show the importance of cooling by 
electron-hole pair creation processes in the presence of strong destructive quantum interfer- 
ence effects and how this cooling mechanism can be used to control these effects by varying 
the temperature of the electrodes. 



II. THEORY 



Model Hamiltonian 



We consider quantum interference effects and vibrationally induced decoherence in elec- 
tron transport through a single molecule that is bound to two metal leads {2 -^]. This 
transport problem can be described by a set of discrete electronic states, which are localized 
on the molecular bridge (M), and a continuum of electronic states, which are localized in 
the left (L) and the right (R) lead, respectively. Thereby, the states on the molecular bridge 



4 



are coupled to the states in the leads. The corresponding model Hamiltonian can be written 



as 



H e l = 2J e mC ] m C m + 2J e kc[c k + (KnfcCfcC m + h.C.) (1) 

m£M fceL.R fceL,R;mGM 

m<ngM 

where the e k denote the energies of the lead states and c k and c k the corresponding creation 
and annihilation operators. Similarly, the mth electronic state on the molecular bridge is 
addressed by the creation and annihilation operators c m and c m , while its energy is given by 
e m . Coulomb interactions may give rise to additional charging energies. These energies are 
accounted for in the Hamiltonian by Hubbard-like electron-electron interaction terms, 
Umn{.c\n c m — <^m) ( c n c n — S n ), where the parameters 5 m distinguish states that are occupied 
(8 m = 1) or unoccupied (8 m = 0) in the molecular system at equilibrium. This distinction 
between occupied and unoccupied states is necessary, because the parameters of the mode l 



Ham iltonian are determined with respect a specific reference state of the junction 



120- 



122j. As a natural choice, we consider the neutral molecule in its electronic ground state 
as a reference system. Thus, for an electronic state above the Fermi level (S m = 0), the 
associated energy e m denotes the energy required to add an electron to the mth electronic 
state of the reference system. For an electronic state below the Fermi level (5 m = 1), 
e m denotes the energy required to remove an electron from state m, in accordance with 
Koopmans' theorem. Note that throughout the paper, the Fermi energy of the leads is set 
to cf = eV. The coupling matrix elements V m k in the third term of H e \ characterize the 
strength of the interaction between the electronic states of the molecular bridge and the 
leads and determine the so-called level-width functions T Kjmn (e) = ^^2 k( z K V^ k V n k5{e — e^) 
(i^=L,R). Modeling the leads as semi- infinite tigh t-bin ding chains with an internal hopping 



parameter 7 = 2eV, these functions are given by 106] 



iWe) = ^^v^T 2 -^-/^) 2 , (2) 



with K G L,R. Thereby, following Refs. 123Nl25j . we assume a symmetric drop of the bias 



voltage $ at the contacts, i. e. the chemical potentials in the left and the right lead are given 
by y^L = e$/2 and /xr, = — e$/2, respectively. 

Applying a bias voltage an electrical current is flowing through the junction. The molec- 
ular bridge may respond to this current, in particular, to the continuous charge fluctuations 



which are induced by the tunneling electrons, by adapting its geometrical structure. As a 
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57|, 



92|-|94, 



126 



resu lt, the vibrational degrees of freedom of the junction may be excited 
128j . The vibrational degrees of freedom are described in our model as harmonic oscillators 
that ar e linearly coup led to the electron (or hole) densities [f^ m c. m — 5 m ) on the molecular 



bridge 
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122 



122| 



(3) 



where the operator a' a denotes the creation operator of the ath oscillator with frequency 
Q a and Q a = a a + a) a the corresponding vibrational displacement operators. The respective 
coupling strengths are denoted by X ma . Note that there is no coupling between the electronic 
states and the normal modes of the molecular junction in the ground-state of the neutral 
molecule, since we employ the normal modes of this state. The role of the parameters 
5 m in the electronic- vibrational coupling term is to impose this constraint. Finally, the 
Hamiltonian of the overall system is given by the sum 



H = H c \ + H vih . 



(4) 



For the nonequilibrium Green's function approach that we use in this article (cf. Sec. Ill Bl) . 
it is expedient to remove the direct electronic-vibrational couplin g te rm i n th e Hamiltonian 
H. To this end, we employ the small polaron transformation 9l|, |92|, |95|, [130] 



H = e s He~ s 



(5) 



E 



kin 



m<n 



with 



S 



A 



X rn = exp[z^-^pp ( 



— fc f c -5 )P 



A, 



(6) 
(7) 



and P a = —i(a a — a^ a ) the momentum operator associated with vibrational mode a. Note 
that, although there is no explicit electronic-vibrational coupling in H, it appears in the 
transformed Hamiltonian H at three different places: i) in the polaron-shifted energies 



e m = e m + (2<5 m — 1) Ea(^L/^a)) H) in additional electron-electron interactions, which add 
to the original electron-electron interaction terms U mn = U mn — 2 J2 a (^ma^na/^a) an d Hi) 
in the molecule-lead coupling term ^ krn iymkX m ^ k c m + h.c), which is renormalized by the 
shift operators X m . 

Because we focus in this article on the effect of electronic-vibrational coupling on the 
transport characteristics of a single-molecule junction, we suppress the effect of electron- 
electron interactions in the following, that is we set U mn = 0. A discussion of the effect of 
electron-electron interactions on the transport characteristics of a molec ular junction can 
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94j . or in Refs. 
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In 



be found, for example, in our earlier work, Refs. 
this context, it is interesting to note that for high bias voltage and within the wide-band 
approximation, electron-electron interactions ha ve no subs tantial influence on the electrical 
current. This was shown by Gurvitz and Prager 



136 



1371 | based on their exact microscopic 



rate equation methodology. At small (finite) bias voltages, however, it has already been 
demonstrated that in the presence of electron-elect ron inter actions quantum interference 
effects become quenched by spin-flip processes [68|, ll38Ml40| . These effects, however, are 
beyond the scope of this work. 



B. Nonequilibrium Green's Function Approach 



To describe charge transport based on the model introduced above, we employ a nonequi- 
librium Green's function approach that was originally proposed by Galperin et al. 9jJ and 

3-941 to account for multiple vibrational modes and 



later extended by Hartle et al. 



51 



multiple electronic states. This scheme allows to describe electronic-vibrational coupling 
nonperturbatively and accounts for the broadening of the electronic states due to the cou- 
pling of the molecular bridge to the leads. It is, therefore, well suited to describe the present 
problem, because quantum interference effects are particularly pronounced in electron trans- 
port through quasi-degenerate electronic states. Moreover, vibrationally induced decoher- 



ence originates from the coupling of such states to the multitude o: 



vibrati onal degrees o f 



freedom of the junction, which, in addition, may be rather strong |21l . l51l . l52l . Il21l Il22l . Il41 |. 

The central quantities in (nonequilibrium) Green's function theory are the single-particle 
Green's functions, which are given for the electronic degrees of freedom by 



G mn {T,r') = -z(T c c m (r)4(r'))^. 



(8) 



They allow to calculate all single-particle observables, such as, e.g., the population of levels 
or the current flowing through the junction. We employ the following ansatz to calculate 
these Green's functions [oi ] : 

G mn (r,r') = -i(T c c m (T)cl(T r )) H (9) 
= -i(T cCm (r)X m (r)4(r')Xt(r')} F 

with the electronic Green's functions G mn (r, r') = — 2(T c c ot (t)c] 1 (t / ))^ and T c the time- 
ordering operator on the Keldysh contour. Thereby, the indices H/ H indicate the Hamilto- 
nian that is used to evaluate the respective expectation values. The effective factorization of 
the single-particle Green's functions G mn into a product of the electronic Green's functions, 
G mn , and a correlation function of shift operators, (T c X m (r)X^(r'))^, is justified, if the dy- 
namics of the electronic and the vibrational degrees of freedom take place at different t ime 



scales. This is conceptually similar to the Born-Oppenheimer approximation 142L Il43| . In 
the present context, however, we do not address the adiabatic but rather the anti-adiabatic 
regime, where the time scales for electron tunneling events are such that the nuclei of the 
molecule can follow the corresponding charge fluctuations. 

Employing the equations of motion for the electronic Green's functions G mn , 

(id T - e m )G mn (r, T')(-id T > - e n ) = 5(r, r')(-id T > - e n ) + S Ljmn (r, r') + S Rjmn (r, r'), 

(10) 

the respective self-energies due to the coupling of the molecule to the left and the right leads, 
S Limn (r, t') and E R mn (r, r'), can be obtained. These self-energies are given to second order 
in the molecule-lead coupling by 

S L /R, m n(r,r')= Yl y^ynk9k{r,r'){T c X n {r')Xl{T)) 1I , (11) 
fceL/R 

where <7fc(r, r') denotes the free Green's function associated with lead state k. The real-time 
projections of these self-energies determine the electronic part of the single-particle Green's 
functions G mn . In the energy-domain the corresponding Dyson-Keldysh equations read 

&) = G^l\e) + G° m t(e) (£ r L ^(e) + S^ p (e)) G^ (e), (12) 

op 

G^(e) = J2 ai mo(t) + ^n,o») °'U^ ( 13 ) 



op 



with 



nrrt n \ 



1 



(14) 



e - e m + irj 

where r\ denotes a positive infinitesimal number. Note that Eqs. (I12p and (|T3l) give the exact 
result in the non-interacting limit, that is for \ ma — > 0, as well as for an isolated molecular 
bridge, where Knfc — > 0. 

Besides the electronic Green's functions G mn , we also need to evaluate the correlation 
functions of the shift operators, (T c X m (r)X^(r'))^, in order to obtain the full single-particle 
Green's functions G mn . To this end, we use a second-order cumulant expansion in the 



dimensionless coupling parameters 
(T c X m (r)Xt(r')) 7r = exp(^z 
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493] 



Vl a Vlp 



D a p{T, T ) - I T^TFl D a p(T, T) ) , (15) 



'aBK 1 > 



with the momentum correlation functions 



= -i(T c P a (r)P^(r')} ff . 



(16) 



Employing the equation of motion for D a p 

' ;-9 T 2 - C^A^r, - fi 2 ) = 5(t, r')6 a p(-d 2 T , - Q 



^20^ 



n, 



el,ay3 



(17) 



we determine the corresponding self-energy matrix U. e \ t ap- The self-energy matrix Tl e \ ia p, 
which describes the interactions between the vibrational modes and the electronic degrees 



of freedom o 
coupling 



91- 



the molecular bridge, is evaluated up to second order in the molecule-lead 



n 



eL,ap{T, T 



Vl a Vlp 



(E mn (r,T')G nm (r',r) + E nm (r' »G ron (r,r')). (18) 



Since H e i, a p depends on the electronic self-energies E mn = £L,mn + ^n,mn and Green's func- 
tions G mn , Eqs. (|TT|) - ([18]) constitute a closed set of coupled nonlinear equations that needs 
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92). 



to be solved iteratively in terms of a self-consistent scheme 

With the Green's functions D a p and G mn , the average vibrational excitation of each 



9 



vibrational mode is obtained according to 921494 ] 



(aU a ) H « -~Im [D< a (t = 0)] - i + £ ^Im[G< m (t = 0)] (19) 
+2 £ ^^Im[G< m (t = 0)]Im[G<„(t = 0)] 

m<n a 

-2 £ ^f^MG< n (t = 0)]lm[G<Jt = 0)], 
where we use the Hartree-Fock factorization 

( C m C mC^ n C n )-ff ~ (c m C m ) (cJjC n ) — (c^C^) (c^ n C m ) (20) 



for m 7^ n. The current is calculated employing the Meir-Wingreen-like formula 



144 1 



92 



7 = 2e / £ E (S^, m „(6)G> m (e) - E> mn ( e )G< m (e)) . (21) 

mn 

III. RESULTS 

In this section we discuss quantum interference and vibrationally induced decoherence in 
single-molecule junctions based on the methodology introduced above. Thereby, we consider 
models of increasing complexity. First, in Sec. IIII Al we discuss basic quantum interference 
effects that occur in electron transport through single molecules. To this end, we investigate 
two complementary model systems, which allows both a comprehensive discussion of the 
most important interference effects and to set up the framework and the methodology that we 
use to analyze quantum interference effects in this article. Next, in Sees. IIII Bl and IIII CI we 
study the effect of electronic-vibrational coupling on the transport characteristics. Thereby, 
we use simplified coupling scenarios, which facilitate the interpretation of the results and 
provide a first understanding of the effect of vibrationally induced decoherence in single- 
molecule junctions. Finally, in Sees. IIII Dl and IIII E\ we consider more realistic electronic- 
vibrational coupling scenarios and discuss the effect of local cooling due to electron-hole 
pair creation processes, which, in the presence of strong destructive quantum interference 
effects, can become dominant (Sec. IIIIDI) . As will be shown in Sec. IIII E|. this effect results 



in a strong temperature dependence of the electrical current, which cannot be explained 



in terms o 



lj I145M148 ) 



the temperature dependence of the Fermi distribution functions in the leads 
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TABLE I: Parameters for the models of single-molecule junctions that are investigated in this 
article. All energy values are given in eV. In all these models the temperature of the electrodes 
Tl+r is set to 10 K. The value of the molecule-lead coupling parameter is v = 0.1 eV. 



The parameters for all models employed are given in Tab. [H The parameters have been 



chosen to represent typical values for molecular junctions simi 
number of first-principles studies of molecular junctions 



21 



48 



ar to those employe d in a 
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122 



141 



1491]. Note 



that we consider only models, where the coupling between the molecule and the electrodes is 
rather weak, that is where the effective factorization of the single-particle Green's functions, 
Eq. (ED, is valid. 



A. Basic Interference Effects 



Quantum interference effects can be constructive or destructive, depending on the model, 
the physical parameters and/or the specific observable. In this section, we study basic 
interference effects for two different models of single-molecule junctions. The first model 
exhibits pronounced destructive interference effects and is referred to as model DES. The 



second model shows constructive interference effects and is ca 
both systems have been extensively studied before 
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led model CON. Note that 



70N82|. They are thus very 



useful to set up the stage for the discussion of vibrationally induced decoherence presented 
in the subsequent sections, Sees. IIII Bl [Til El It is also noted that similar models have been 
used to study interference effects in quantum dot arrays, in particular in the context of 



Aharonov-Bohm interferometers 



69 



150 



152 1 . 



Interference effects arise in electron transport through molecular junctions whenever elec- 
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FIG. 1: Panel (a): Schematic representation of example wavefunctions of the electronic eigenstates 
of model DES. Panel (b): Localized molecular orbitals corresponding to the eigenstates of model 
DES. Panel (c): Schematic representation of example wavefunctions of the electronic eigenstates 
of model CON. Panel (d): Localized molecular orbitals corresponding to the eigenstates of model 
CON. 

trons can traverse the junction along different but equivalent pathways. This is the case, for 
example, if the electrical current of a junction is carried by two quasidegenerate electronic 
states. The first model system that we study, model DES (see Tab. |T] for a summary of the 
respective parameters), comprises two electronic states, which are located at e± = 0.5 eV and 
€2 = 0.505 eV above the Fermi level of the junction. Both of these states are coupled to the 
left lead with coupling strengths — 0.1 eV and to the right lead with coupling strengths 
^r.i = 0.1 eV and z/ Rj 2 = — 0.1 eV, respectively. These coupling strengths, in particular the 
different sign in the coupling to the right lead, reflect the different spatial symmetry of the 
two states, which represent symmetric and antisymmetric combinations of localized molec- 
ular orbitals (see Figs. [Qi and [TJd for a graphical representation of the two electronic states 
and the corresponding localized molecular orbitals). The coupling to the electrodes induces 
a broadening of the two states of ~ 20 meV that exceeds their level spacing — £i = 5 meV, 
i.e. they are quasidegenerate. 

The current- volt age characteristic of junction DES is represented in Fig. [3] by the solid 
purple line. Due to destructive quantum interference effects, it shows rather low levels of 
electrical current. This is demonstrated by comparison with the dashed purple line, which 
depicts the current-voltage characteristic of this junction disregarding quantum interference 
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(a) (b) (c) (d) 



LMR LMR LMR LMR 

e 2 




FIG. 2: Panel (a) and (b): Resonant and non-resonant electron transport processes through junc- 
tion DES. The purple wiggly lines depict the coherent sum of the outgoing wavefunctions that 
are associated with electron tunneling through each of the two electronic states. Panel (c) and 
(d): Resonant electron transport through state 1 and 2 if the other state (i.e. state 2 and 1, re- 
spectively) would not be present. Due to destructive interference, the corresponding tunneling 
amplitudes (purple wiggly lines) are much larger than their coherent sum (depicted in Panel (a)). 

effects, that is, the incoherent sum of the electrical currents that are flowing through either of 
the two electronic states. As can be seen, quantum interference effects suppress the current 
level in this junction by more than an order of magnitude. This applies throughout the 
whole range of bias voltages considered, including the non-resonant as well as the resonant 
transport regime, which are separated by the step that appears at e$ ~ 2ex/2 in both 
characteristics. 

Quantum interference effects are destructive in this system, because the outgoing wave- 
function, which is associated with an electron tunneling event through state 2, differs by a 
phase of 7r from the the one, which is associated with tunneling through state 1. This is 
a result of the different spatial structure of the two states and is reflected in the different 
signs of the respective molecule-lead coupling strengths to the right lead: z/ R1 = — u R 2 - 
The two outgoing wave-functions associated with electron tunneling through state 1 and 2, 
therefore, destructively interfere with each other, leading to a strong suppression of the re- 
spective tunnel current. This is illustrated by an example for a resonant and a non-resonant 
transport process in Figs. Ek and Eb, respectively, where the coherent sum of the outgoing 
wavefunctions is depicted. Figs. Eb and Eli show the two outgoing wave-functions associated 



13 



5 
4 

< 
^3 



a 

— 

o 



■ with interference 
w/o interference 
■|v/.//U/2l=0.07eV 



: 







0.5 1. 1.5 

bias voltage $ [V] 



FIG. 3: (Color online) Current-voltage characteristics of the linear molecular conductor described 
by model DES (see Figs.QJt and[tb). The solid purple line depicts the current- volt age characteristic 
of this junction, including quantum interference effects. The dashed purple line is obtained by 
discarding them. If the absolute value of the coupling strengths z^l/r,i/2 of model DES is reduced 
to 0.07 meV, the current- voltage characteristic depicted by the solid turquoise line is obtained. 

with the resonant transport process in Fig. [2^. Note that these interference effects are more 
pronounced the closer the two electronic states are in energy. 

Employing a different representation, the suppression of the current level in model DES 
can also be understood on different grounds. Thereby, the two eigenstates of this system 
are unitarily transformed to two localized states (see Figs. [T^i and [lb for a schematic repre- 
sentation of the two equivalent representations). These states have the same energy eo and 
are coupled with each other by the coupling strength A = (e 2 — ei)/2. They represent, for 
example, the left and the right part of a molecular conductor. Using this local represen- 
tation, electron transport through this junction involves a sequence of processes: electron 
transfer from the left electrode to left part of the molecule, followed by an intramolecular 
electron transfer process from the left part of the molecule to the right part, and finally a 
transfer process from the molecule to the right electrode. For | z y L/R,i/2 1 ^ the bottleneck 
for electron transport through this linear molecular conductor is the intramolecular electron 
transfer process. Accordingly, the suppression of the electrical current flowing through junc- 
tion DES can also be understood in the local representation by the small effective coupling 
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A rather than in terms of quantum interference effects. Other aspects of the problem are, 
however, more straightforwardly understood in the eigenstate representation. For example, 
for r L ,iir Rjll > A 2 , a reduction in the coupling to the leads results in an increase of the 
current that is flowing through junction DES. This is illustrated by the solid turquoise line 
in Fig. [21 which shows the current-voltage characteristic of model DES with a reduced cou- 
pling to the leads, ^l,i/2 = vr,i = — ^r,2 = 0.07eV. While this is easily understood in terms 
of a suppression of destructive interference effects due to a less pronounced broadening of 
the electronic states, it may hardly be anticipated using the localized state picture, where 
transport is understood rather as a sequence of intramolecular electron transfer processes. 

To further analyze quantum interference effects in model DES and facilitate the discussion 
of more complex models, we e mploy the transmission function t(e) of model DES, which is, 



in general, defined by 



144 1 



t(e) = Yl r L , mn (e)G^(6)r R>op (e)G^(e) (22) 

mnop&A 

= 4vr 2 \tkk>(e)\ 2 5(e-e k )5(e-e kl ). (23) 



fcgR.fc'eL 



Thereby, 



M<0 = E^mnKfe* (24) 



denote the transition matrix elements that represent the transmission amplitudes of the 
conduction process. Accordingly, the transmission function can be split into an incoherent 
term 

t inc (e)=4vr 2 £ J2KkG^ n V nk/ \ 2 5(e-e k )5(e-e k/ ), (25) 

fceR,fc'eL mn 

^2 r2 



\e — ei + zlp \e — £ 2 + iTf 
which represents the incoherent sum of the transmission amplitudes V^ k G T mn V n k' , and an 
interference term 

Me) = *(e) -Me)> ( 27 ) 
= -2Re — ^ — , (28) 

which encodes the respective interference effects. Thereby, to derive the expressions ( |26|) 
and f f28|) . we have employed the wide-band approximation Y = T Lt n(fi L ). The structure of 
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FIG. 4: (Color online) Panel (a): The solid black line shows the transmission function of the 
linear molecular conductor DES (graphically illustrated in Figs. [T^ and [lb)- In addition, the 
corresponding incoherent (dashed orange line) and interference term (solid green line) are depicted. 
Panel (b): The same as Panel (a), where, however, the energy of the second electronic state is shifted 
from e 2 = 0.505 eV to e 2 = 0.53 eV. 

these expressions suggests an equivalent decomposition of the transmission function t(e) in 
terms of transmission amplitudes Au 2 



t(e) = |A!-A 2 | 2 , 
t inc (e) = |A!| 2 + |A 2 | 2 , 
t int (e) = A*A 2 + A^Ai, 



(29) 
(30) 
(31) 
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with 

Ai = (32a) 

e — Ei + tL 

A 2 = ] —. (32b) 

e - e 2 + iT 

These amplitudes describe electron transport through state 1/2 if the other state 2/1 is not 
present, i.e. in the limit e 2 /i — > oo. It should be noted that the definition of the incoherent 
term fl25J) and the interference term f|2T|) is not unique, and depends, for example, on the basis 
that is used to represent the electronic degrees of freedom of the molecular bridge. Other 
decomposition schemes, which employ, e.g., molecular conductance orbitals, are possible 
and have been used to study interference effects 74j . In this work, however, we seek for an 
understanding of the transport properties of a molecular conductor in terms of its eigenstates, 
which, in contrast to molecular conductance orbitals, represent an energy-independent basis. 
For later reference, it is noted that the decomposition of the transmission function t(e) into 
the incoherent term ( 125]) and the interference term fT271) can also be obtained by disregarding 
the off-diagonal elements of the self-energy matrix ' y in the current formula fl2T]) (cf. App. 

E). 

Fig. shows the transmission function of junction DES (solid black line) as well as 
the corresponding incoherent (dashed orange line) and interference term (solid green line). 
Due to the quasi- degeneracy of the two states, the interference term is almost as large as 
the incoherent term. It shows only negative values, which indicate destructive interference 
effects. Thus, the incoherent term, which describes electron transport through this system 
without quantum interference effects, and the interference term almost cancel each other, 
leaving only a small peak at e w (e 2 — ei)/2 in the respective transmission function. This 
peak is associated with the low current levels that are obtained for this model system (cf. 
Fig. E]). 

Increasing the level spacing in this model by shifting the energy of the second state to 
higher values, e 2 = 0.53 eV, the tranmission function becomes more complex, as shown in 
Fig. HI The interference term is no longer restricted to negative values, but turns from 
negative to positive values if the energy of the electron is located between the energy of 
the first and the second state. Hence, for a larger level spacing, junction DES exhibits 
both destructive and constructive interference effects. This result is more straightforward 
understood based on the interference picture in the eigenstate representation than in terms 
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FIG. 5: (Color online) Current-voltage characteristics of the branched molecular conductor de- 
scribed by model CON (see Figs, [lb andQJl). The solid purple line depicts the current-voltage 
characteristic of this junction, including quantum interference effects. The dashed purple line is 
obtained by discarding them. 

of the intramolecular coupling A in the local representation, especially if these constructive 
interference effects are dominant, such as, for example, if e\ < ep < e 2 . 

Next, we consider model CON. It differs from the previous system only by the sign of the 
coupling strength z/ R 2 (see Tab.flJ. The two electronic states of model CON are thus coupled 
to the right lead by the same coupling strengths (z/ R1 = z/ R 2 )- Due to the different sign, 
model CON and model DES describe very different types of molecular conductors. While 
model DES corresponds in the local representation to a linear molecular conductor, model 
CON describes a branched molecular conductor (cf. Figs. [Tfc andCQi). Similar to model DES, 
it can also be unitarily transformed and represented by two localized molecular states that 
are mutually coupled with each other with the coupling strength A = (e 2 — ei)/2. While 
in model DES, however, each of these states is coupled to a different lead, model CON 
can be mapped onto a delocalized state, which corresponds to the backbone of a molecular 
conductor that is connected to both electrodes, and a localized state, which corresponds to 
a side-branch of the molecular conductor that is not directly connected to the electrodes of 
the junction. 

The current-voltage characteristic of this junction is represented in Fig. [5] by the solid 
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FIG. 6: (Color online) The solid black line shows the transmission function of the branched molec- 
ular conductor CON (graphically illustrated in Figs, [lb andHH). In addition, the corresponding 
incoherent (dashed orange line) and interference term (solid green line) are depicted. 

purple line. In contrast to model DES, it does not show pronounced quantum interference 
effects. The current flowing through this junction is almost the same as the incoherent 
sum of the currents that are flowing through each of the two states. This can be inferred 
by comparison with the dashed purple line, which depicts the respective incoherent sum 
current. The only difference between these two characteristics is an enhanced broadening of 
the step that separates the non- resonant and the resonant transport regime at e$ ~ 2ei/ 2 . 

To understand this result, we analyze the corresponding transmission function. As for 
model DES (cf. Eqs. 1221) . it can also be decomposed into an incoherent and an interference 
term (cf. App. IB"I) : 

i inc (e) = |A 1 (e)| 2 + |A 2 (e)| 2 , (33) 
t int (e) = A*(e)A 2 (e)+A*(e)A 1 (e). 

with the transmission amplitudes 

Ai = ■ T ^ 2M1 _ e2 , (34) 

e — ei + iT 1 2 

A 2 = , ^EEEE ' ( 35 ) 

e — e 2 + u 1 2 
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which, similar to the amplitudes ( 132|) . correctly describe electron transport through junction 
CON in the two limits €2/1 — > 00. Fig. [6] depicts the transmission function of this system 
and the corresponding incoherent and interference term. For low energies e < €\ as well as 
high energies e > €2, constructive interference effects, which are associated with a positive 
interference term, are indeed active and double the transmittance of the junction. This 
leads to significantly larger current levels in the non-resonant transport regime, i.e. for 
e$ < 2ei (cf. Fig. [5]). For intermediate energies e\ < e < 62, however, the interference term 
indicates a transition from constructive to destructive interference effects, similar as for an 
increased level spacing in model DES (see Fig. Hb). This is accompanied by a steep increase 
and decrease of the incoherent and the interference term, respectively. At e = (t2 + ei)/2, 
the two terms are exactly the same and cancel each other. The transmission function of 
model CON thus drops to zero at this energy, exhibiting an antiresonance. These strong 
destructive interference effects reduce the increase of the current level at the onset of the 
resonant transport regime (e$ ~ (e2+ei)), leading, in conjunction with the increased current 
level in the non-resonant transport regime, to a much broader step in the current-voltage 
characteristic. At higher bias voltages, destructive and constructive interference effects 
cancel each other such that the same current level is obtained as if interference effects would 
play no role in this system. 

B. Vibrationally Induced Decoherence in the Resonant Transport Regime 

In this section, we study the influence of electronic-vibrational coupling on the quantum 
interference effects that we found in the resonant transport regime of model DES and model 
CON {i.e. for bias voltages e$ > 2e!/ 2 , cf. Sec. IIII Al) . The influence of this coupling on the 
corresponding effects in the non-resonant transport regime (e$ < 2ei/2) are addressed in 
Sec. lfflCl 

To study vibrational effects in model DES and model CON, we extend these model 
systems by coupling to a single vibrational degree of freedom. Thereby, we consider first 
scenarios, where the vibrational mode is coupled to one of the states only, An = OeV 
and A21 = 0.05 eV. The resulting model systems are referred to as model DESVIB and 
model CONVIB (see Tab. [J for the complete set of parameters). The specific choice of the 
electronic-vibrational coupling strengths in these systems is used to simplify the discussion of 
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the corresponding decoherence mechanisms, where, in addition, the energy of the electronic 
states in models DESVIB and CONVIB have been chosen such that the polaron shifted 
levels e m agree with those of models DES and CON. This allows to separate static effects 
of the vibrations, due to a polaron shift of the energies, from dynamical decoherence effects 
(vide infra). Systems with An 7^ A21 show very similar decoherence effects. Such s yst ems 
are discussed, for example, in Ref. |52j and in Sees. IIIIDI and lHIEl Note that in Ref. 52j we 
have already studied and analyzed the transport characteristics of a single- molecule junction, 
which is very similar to junction DESVIB. 

The current- voltage characteristic of model DESVIB is represented in Fig. [7|l by the solid 
black line. It is referred to as the vibronic current-voltage characteristic in the following. 
The corresponding electronic current- voltage characteristic, where the electronic-vibrational 
coupling of model DESVIB is discarded, is depicted by the solid purple line. Both current- 
voltage characteristics show a step at e$ = 2ex/ 2 and e$ = 2ei/ 2 , respectively, which indicates 
the onset of resonant transport processes (cf. Sec. IIII Al) . Due to the polaron-shift of the 
two electronic states, this occurs at slightly lower bias voltages in the vibronic current- 
voltage characteristic than in the electronic current- volt age characteristic. In contrast to 
the electronic case, however, the vibronic current- voltage characteristic exhibits a number of 
additional steps at higher bias voltages, e$ — 2( e-|/o+nfM (n EN). These steps correspond 
to the onset of resonant excitation processes 57|, |92J, where electrons tunnel resonantly onto 



the molecular bridge, exciting the vibrational mode by n vibrational quanta. An example 
of such an excitation process is graphically depicted in Fig. Ek- 

The comparison of the results in Fig. [TJ in particular with respect to the current- voltage 
characteristic of model DES (dashed purple line), show that electron- vibrational coupling 
results in a pronounced increase of the current level. This is caused by vibrationally induced 
decoherence that originates from the different vibronic coupling strengths of the two elec- 
tronic states. For this particular junction, where An = eV, electrons traversing the junction 
can only interact with the vibrational degree of freedom if they pass through state 2. The 
interaction with the vibrational mode in this state thus provides 'which-path-information' 



that destroys the coherence of the electrons 



52]. 



Without these decoherence effects, the current levels of model DESVIB would be an order 
of magnitude smaller and, actually, be the same as for model DES. This is not straightfor- 
wardly evident from the corresponding model parameters, since the level spacing in model 
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(a) 





FIG. 7: (Color online) Panel (a): Current- voltage characteristics of the linear molecular conductor 
DESVIB (solid black line). The solid purple line shows the current-voltage characteristic of this 
junction disregarding electronic-vibrational coupling. The dashed purple line shows the current- 
voltage characteristic of model DES. The dashed gray and red lines depict the current-voltage 
characteristic of junction DESVIB that is obtained with the vibrational degree of freedom kept in 
thermal equilibrium at 10 and 5000 K, respectively. Panel (b): Vibrational excitation characteristic 
corresponding to the vibronic current-voltage characteristic shown in Panel (a). 
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FIG. 8: Examples for resonant inelastic electron transport processes, where the vibrational degree 
of freedom of the junction is singly excited (Panel (a) and (b)) and deexcited (Panel (c) and (d)) 
upon an electron transfer process from the left electrode onto the molecular bridge (Panel (a) and 
(c)) or from the molecular bridge to the right lead (Panel (b) and (d)). 

DESVIB (e 2 — t\ = 30meV) is larger than in model DES (e 2 — e% = 5meV). Accordingly, 
the electronic current- voltage characteristic of model DESVIB (solid purple line in Fig. [7^) 
exhibits much higher current levels than the one of model DES (dashed purple line in Fig. 
EK, which is the same as the solid purple line in Fig. [3]). The vibronic current- voltage char- 
acteristic (solid black line in Fig. [7|l) displays even larger current values. This, however, 
cannot be attributed to a different level spacing, since the polaron-shifted energy levels Z\ 
and €2 of model DESVIB coincide with the energy levels of junction DES. Thus, on the 
level of the electronic parameters, destructive quantum interference effects should be as pro- 
nounced in model DESVIB as in model DES. As this is not the case and the current level of 
model DESVIB even reaches the same values as the incoherent sum current of model DES 
(dashed purple line in Fig. [3]) at higher bias voltages ($ > 1.5 V), we conclude that in model 
DESVIB electronic- vibrational coupling strongly quenches destructive quantum interference 
effects and that this quenching becomes almost complete at higher bias voltages. 

The significantly enhanced electrical current of junction DESVIB is not the only mani- 
festation of vibrationally induced decoherence in this system. Also, the relative step heights 
in the corresponding current-voltage characteristic are strongly influenced by quantum in- 
terference effects and vibrationally induced decoherence. Typically, the relative step heights 
are correlated with the Franck-Condon matrix elements Fq u = exp(— A^/f^XAji/^i™)/ 71 " 
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These matrix elements determine the probability for a transition from the vibrational ground 
state to the nth excited state, which, for example, may be part of an inelastic excitation 
process (Figs. Ek). Indeed, there is no strict one to one correspondence between the step 
heights and the Franck-Condon factors F 0n , because other processes, like the ones depicted 
in Figs. [Hb - Fig. [SH, and the respective nonequilibrium excitation of the vibrational mode 



lead, in general, to a suppression of the current level |57j. But the first steps in the vibronic 
current-voltage characteristic of model DESVIB are, nevertheless, far too low. Instead of a 
ratio of 1 : 1.2 between the first and the second step, as it is expected for transport through 
a single level {5?]], one observes a ratio of 1 : 2.2. The successive steps show a similar behav- 
ior. This indicates that at the onset of the resonant transport regime destructive quantum 
interference effects decrease the current level of this system and that, at higher bias voltages, 
these effects become gradually suppressed as more inelastic processes become active. This is 
due to an increase of both the number of active excitation and deexcitation processes (Figs. 

- EH), where the latter is particularly enhanced by the corresponding level of vibrational 
excitation (depicted in Fig. [7b ) that increases rapidly with the applied bias voltage in the 
resonant transport regime. 

As before in Sec. IIII A\ these findings can be analyzed in more detail, considering the 
transmission function of this system. To define a transmission function in the presence of 
electronic-vibrational coupling, we restrict the vibrational degree of freedom to its thermal 
equilibrium state at temperature T, employ the wide-band approximation and consider high 
bias voltages (e$ ^> ei/2). With these assumptions, the current through junction DESVIB 
can be expressed in a form similar to Landauer theory 
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I^2e (de/2vr)t(e), (36) 

J —00 

with the transmission function 

t(e) = ztr{r L G > }. (37) 

A decomposition of this transmission function in terms of transmission amplitudes is not 
obvious. Instead, we use the observation (cf. Sec. IIII Ap that the incoherent term of the 
transmission function can also be obtained by neglecting the off-diagonal elements of the 
self-energy matrix in the current formula (I2T]) . The interference term is then given by 
£mt( e ) = —tmc( e )- For our specific model system, we thus obtain (see App. |A]for details) 
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FIG. 9: (Color online) The solid black lines show the transmission function of the linear molecu- 
lar conductor DESVIB for an effective temperature of the vibrational degree of freedom of 10 K 
(where no offset is added) and 5000 K (where an offset of 2 is added to these lines). In addition, 
the corresponding incoherent (dashed orange lines) and interference terms (solid green lines) are 
depicted. 



oo -p 2 

2 + A 2 £ Z,(s)eW _ 1 (38a) 



|e — ei + iT\ , \e — e 2 — + iT 

I I i = — oo 1 



t int (e) = -2A 2 r 2 Re 



1 



(38b) 



(e-e 2 + ir)(e-ei-ir). 

Thereby, the prefactor A = e _ ^2i/( 2n i))( 2Ar vib+i) j s determined by the average vibrational 
excitation iV v ib = (e^ 1 — 1) _1 and the inverse temperature (3 = (^bT 1 ) -1 , while h{x) = 



(2(A| 1 /S7f ) a/ N vih (N vih + 1)) denotes the Zth modified Bessel function of the first kind. 
The transmission function of model DESVIB obtained in this way is depicted in Fig. [9] by 
the solid black line. The corresponding incoherent and interference terms are given by the 
dashed orange and solid green line, respectively. Thereby, two scenarios are distinguished: 
(i) the effective temperature of the vibrational degree of freedom is assumed to be 10 K, i.e. 
it is effectively not excited, and (ii) the temperature is assumed to be 5000 K, which results 
in a level of vibrational excitation that is comparable to the one obtained in the resonant 
transport regime (cf. Fig. [7b). In both cases, the interference term is significantly weaker 
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with respect to the incoherent term than in the respective electronic transport scenario, 
which, due to the polaron-shift of the electronic energy levels, is the one of model DES 
(cf. Fig. Hk. for the respective transmission function and the corresponding incoherent and 
interference term). Accordingly, the transmittance of the junction is also larger and the 
corresponding current levels are higher. This can bee seen in Fig. [T^l by comparison of the 
dashed gray and red line, which depict the current- volt age characteristic of model DESVIB 
with the vibrational mode kept in thermal equilibrium at 10 and 5000 K, respectively, with 
the dashed purple line, which depicts the electronic current- volt age characteristic, including 
the polaron shift of the electronic energy levels. 

The suppression of the interference term indicates a strong quenching of quantum inter- 
ference effects in this system. This quenching originates from both a static mechanism that 
constitutes an effective renormalization of the level-width functions r# jTOm and vibrationally 
induced decoherence. The factor A 2 , which precedes the interference term (138bl) . includes 
both of these effects. It describes, on one hand, the suppression of electronic tunneling events 
(as, e.g., depicted in Fig. [2]) due to electronic- vibrational coupling or, equivalently, the de- 
creased overlap between the vibrational states of different charge states of the molecular 
bridge. At 10 K, where the vibrational mode is effectively restricted to its ground state, this 
suppression corresponds to the Franck-Condon matrix element Fqq. On the other hand, the 
prefactor A 2 also includes the effect of vibrationally induced decoherence, which originates 
from inelastic tunneling processes, where electrons tunnel resonantly through the junction 
at energies ~ ei/2- Examples for such processes are depicted in Figs. [8b and EH. Inelastic 
processes, where electrons tunnel resonantly through the junction at energies e ~ €2 + nfli 
with n G Z/{0} (see Figs. [Hh^ and [Hb), complement the effect of vibrationally induced deco- 
herence and result in additional side-peaks in the transmission function and the incoherent 
term but not in the interference term. Both the suppression of the interference term by the 
prefactor A 2 and the vibrational side-peaks become stronger at higher levels of vibrational 
excitation, which increases the probability for inelastic processes, in particular deexcitation 
processes. As a result, the main peak in the interference term, which appears at e ~ ex/2 
and is clearly visible in the electronic case, is somewhat reduced in the vibronic case with 
10 K and vanishes almost completely in the vibronic case with 5000 K. Accordingly, side- 
peaks appear only in the vibronic scenarios and are more pronounced for a higher effective 
temperature of the vibrational mode. Thus, vibrationally induced decoherence leads to a 
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complete suppression of interference effects in this system as the effective temperature of 
the vibrational mode increases. 

Next, we consider model CONVIB. The corresponding current- volt age characteristic and 
the respective level of vibrational excitation are shown in Fig. (TDJ In contrast to model 
DESVIB, these characteristics do not exhibit pronounced decoherence effects, because in 
model CON interference effects play, a priori, a less significant role in the resonant transport 
regime (cf. Sec. Ill I Aj) . Moreover, the vibrational excitation characteristics of model CON- 



VIB and model DESVIB almost coincide, since inelastic processes, due to the simplified 
electronic-vibrational coupling scenario with An = 0, are not influenced by interference ef- 
fects in both systems. Analysis of the corresponding transmission function, however, reveals 
that electronic-vibrational coupling affects quantum interference effects that are active in 
model CONVIB in a very similar way as in model DESVIB. 

Using the same scheme as before, the transmission function and the interference term of 
model CONVIB are given by (see App. [B| 

t( e ) = c(e) [\e-e 2 + iT\ 2 + 2A 2 {e-e 1 ){e-e 2 )-2A 2 T 2 (2-A 2 )] (39a) 

oo 

+A 2 h(x)e^ 2 c(e - m x ) (|e -e x - m x +iV\ 2 - T 2 A 2 ) , 

l=—oo 

t iQt (e) = c(e)[2A 2 (e-e 1 )(e-e 2 )-2A 2 T 2 (l-A 2 )], (39b) 



with 



c e 



(40) 



(e-e! + ir)(e-e 2 + zr) + A 2 r 2 

They are represented in Fig. [Ill including the respective incoherent term. Vibrationally 
induced decoherence causes the same effects as outlined before, that is, a suppression of 
the interference term and the appearance of side-peaks that are not counterbalanced by the 
interference term. As a result, the antiresonance at e = (?i +?2)/2 is quenched and even 
vanishes as the level of vibrational excitation increases. 



C. Vibrationally Induced Decoherence in the Non-Resonant Transport Regime 

While in Sec. IIIIBI we have focused on the resonant transport regime of a molecular 
junction, here we investigate the effect of electronic-vibrational coupling on quantum in- 
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FIG. 10: (Color online) Panel (a): Current- voltage characteristics of the branched molecular con- 
ductor CONVIB (solid black line). The solid purple line shows the current-voltage characteristic 
of this junction disregarding electronic- vibrational coupling. The dashed purple line shows the 
current-voltage characteristic of model CON. The dashed gray and red lines depict the current- 
voltage characteristic of junction CONVIB that is obtained with the vibrational degree of freedom 
kept in thermal equilibrium at 10 and 5000 K, respectively. Panel (b): Vibrational excitation 
characteristic corresponding to the vibronic current- volt age characteristic shown in Panel (a). 
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FIG. 11: (Color online) The solid black lines show the transmission function of the branched 
molecular conductor CONVIB for an effective temperature of the vibrational degree of freedom of 
10 K (where no offset is added) and 5000 K (where a horizontal offset of 1.5 and a vertical offset of 
0.05 is added to these lines). In addition, the corresponding incoherent (dashed orange lines) and 
interference terms (solid green lines) are depicted. 

terference effects in the non-resonant transport regime. To this end, we discuss the trans- 
port characteristics of the same model systems as in Sec. IIIIBI but for lower bias voltages: 
e$ < 0.25 eV< 2e 1/2 . 

The conductance- voltage characteristic (i.e. the differential conductance d//d$) of model 
DESVIB is depicted in Fig. H2] by the solid black line. For very small bias voltages, e$ < fli, 
the conductance of junction DESVIB is almost constant. This is related to the fact that 
electrons can tunnel through the junction only in non-resonant, electronic tunneling events 
(as depicted in Fig. |2b). As the number of these transport processes increases linearly 
with the applied bias voltage, and because these processes are almost equally probable, the 
associated current level increases also almost linearly, giving rise to the almost constant level 
of conductance. At higher bias voltages e$ > Qi, the conductance of the junction increases 
substantially. Here, electrons have enough energy to excite the vibrational degrees of freedom 
upon transport through the junction. An example of such a process is depicted in Fig. [T3k . 
These processes enhance the transmission probability for electrons that have enough energy, 
that is an energy in the range [/xr+^i,/^l] (for e$ > Q\). Accordingly, the conductance of the 
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FIG. 12: (Color online) Conductance-voltage characteristics of the linear molecular conductor 
DESVIB. The low-bias regime (<£ < 0.1 V) is highlighted in the inset of this figure, where, in 
addition, a logarithmic scale is employed. Thereby, the solid purple line shows the current- volt age 
characteristic of this junction disregarding electronic- vibrational coupling. The dashed purple line 
shows the current-voltage characteristic of model DES. The dashed gray line depicts the current- 
voltage characteristic of junction DESVIB with the vibrational degree of freedom kept in thermal 
equilibrium at 10 K. 

junction increases, leading to a step in the conductance- voltage characteristic at e$ > Oi . 
Such steps occur whenever the bias voltage matches a multiple of the vibrational frequency, 
e$ ±2 nQi (n G N), corresponding to the opening of a new inelastic transport channel, 
where the vibrational degree of freedom is excited by n vibrational quanta. This behavior 
is well known and used in inelastic electron tunneling spectroscopy of molecular junctions 



m 



47H50J. Note that, despite the opening of inelastic excitation channels at e$ — nQ± 



(n £ N), the corresponding level of vibrational excitation is rather low in this regime. This 
is due to non-resonant (cf. Figs. [T3b and [T3b ) as well as resonant deexcitation processes (cf . 

94| . which cool the vibrational degree of freedom in this bias regime 



Figs. Efc, Eli and M) 
very effectively. 

The effect of vibrationally induced decoherence in this regime is, similar to the resonant 
transport regime (cf. Sec. IIIIBl) . a pronounced increase of the current. This is demonstrated 
in Fig. [12] by comparison of four different conductance- voltage characteristics: the vibronic 
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FIG. 13: Examples of non-resonant inelastic processes, including transport (Panel (a) and (b)) and 
electron-hole pair creation processes (Panel (c)). Note that transport processes can occur at low 
temperatures (i.e. /cb^l+r ^> ^1) via the excitation (Panel (a), e<3> > f^i) or deexcitation of the 
vibrational degree of freedom (Panel (b)) while, in contrast, electron-hole pair creation processes 
involve predominantly deexcitation processes, because the respective excitation processes would be 
suppressed by Pauli blocking. 

(solid black line) and the electronic conductance- volt age characteristic of model DESVIB 
(solid purple line), the electronic conductance- voltage characteristic of model DES (dashed 
purple line) and the conductance-voltage characteristic of model DESVIB with the vibra- 
tional degree of freedom kept in thermal equilibrium at 10 K (dashed gray line). Overall, 
we find the same three mechanisms at work as in the resonant transport regime (cf. Sec. 
IIIIBI) . including vibrationally induced decoherence due to inelastic processes. First, the 
comparison of the electronic conductance-voltage characteristics of junction DESVIB and 
DES shows that destructive quantum interference effects are enhanced by the polaron-shift 
of the two electronic states due to an increase of their quasidegeneracy (note that the energies 
of the electronic states in model DES correspond to the polaron-shifted energies in model 
DESVIB). Second, the suppression of electronic tunneling events through state 2, which 
is due to the coupling of this state to the vibrational mode and which is, in this regime, 
approximately given by the Franck-Condon factor Fqq, leads to a significant quenching of 
destructive interference effects. This can be seen by comparison of the conductance levels 
for the vibronic and the electronic cases, in particular at bias voltages e$ < where only 
electronic tunneling events occur. 
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Besides these static effects, for voltages above the first inelastic threshold e$ > fii, there 
is a third dynamical decoherence effect due to non-resonant excitation and deexcitation 
processes (see Figs. [T3"a and [T3"b). These processes can only occur via (virtual) tunneling 
processes through state 2, because only this state is coupled to the vibrational degree of 
freedom. Hence, these processes are not influenced by the interference effects in this system. 
In contrast, non-resonant electronic transport processes (Fig. [2b) are suppressed by destruc- 
tive interference effects in this system and, consequently, their contribution to the current is 
rather small. This can be inferred by comparing the conductance values at e<3> < Q\, where 
only electronic transport processes contribute, and flj < e$ < 2f2 1; where inelastic processes 
become active. The ratio of these values is approximately 1 : 10. Calculations, where we 
discard the destructive interference effects in this system, give a much smaller ratio of these 
conductance values, that is, approximately 1 : 1.3. Therefore, non-resonant excitation and 
deexcitation processes increase the current level of the junction in the non-resonant trans- 
port regime in a very similar way as resonant excitation and deexcitation processes do in 
the resonant transport regime (cf. Sec. IIIIBI) . 

As was shown in Sec. lIIIBl vibrationally induced decoherence is thereby more pronounced 
the higher the vibrational degree of freedom is excited. This can be seen by comparison of 
the vibronic conductance- voltage characteristic (solid black line), where the finite level of 
excitation of the vibrational mode is taken into account, with the conductance- voltage char- 
acteristic that is obtained by restricting the vibrational degree of freedom to its thermal 
equilibrium state at 10 K, that is, effectively, to its ground state (dashed gray line). Overall, 
the higher level of vibrational excitation in the vibronic case leads to significantly larger con- 
ductance values. Moreover, as the number of inelastic processes and, thus, the corresponding 
level of vibrational excitation increases almost linearly with the applied bias voltage, the con- 
ductance of the junction also increases almost linearly above the inelastic threshold e$ > Q±. 
Note that the vibronic conductance-voltage characteristic displays higher conductance val- 
ues even for e$ < Q±. This is due to higher order effects in the molecule- lead coupling that 
result not only in a finite population of the electronic level but, in the same way, also in a 
finite level of vibrational excitation ((a^a) « 0.05) at these bias voltages. 

Next, we investigate non-resonant transport through model CONVIB, where we con- 
sider the same vibronic and electronic transport scenarios as before. The correspond- 
ing conductance-voltage characteristics are depicted in Fig. [TJ] for small bias voltages 
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FIG. 14: (Color online) Conductance-voltage characteristics of the branched molecular conductor 
CONVIB. The solid purple line shows the current-voltage characteristic of this junction disregard- 
ing electronic-vibrational coupling. The dashed purple line shows the current- voltage characteristic 
of model CON. The dashed gray line depicts the current-voltage characteristic of junction CONVIB 
that is obtained by evaluating the vibrational degree of freedom in the thermal equilibrium state 
it would acquire at 10. 

($ < 0.25 V). As in model DESVIB, the polaron-shift of the electronic levels enhances quan- 
tum interference effects. The conductance associated with electronic transport in model 
CON is therefore larger than in model CONVIB. In contrast, the vibronic conductance- 
voltage characteristic exhibits overall smaller conductance values. This is due to both a 
reduction of the probability for electronic transfer events (Fig. [2b) as well as vibrationally 
induced decoherence, which quenches the constructive interference effects that are active in 
this system. Thereby, a higher level of vibrational excitation, due to a more important role 
of deexcitation processes, leads to a further reduction of conductance (compare the solid 
black and the dashed gray line). The steps, which appear very pronounced in the vibronic 
conductance- voltage characteristic of model DESVIB at e$ ±2 nfl% (n G N), are, therefore, 
substantially reduced, which reflects the fact that at these bias voltages the reduction of 
conductance due to vibrationally induced decoherence and the increase of conductance due 
to the onset of non-resonant excitation processes compete with each other. Another effect of 
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vibrationally induced decoherence is that the slope of the vibronic conductance- voltage char- 
acteristic can even be negative after the onset of an inelastic channel at e$ ~ nQi (n EN). 
This can be seen, for example, after the first step in the solid black line of Fig. [T3] and is 
related to the number of inelastic processes, which increases with the applied bias voltage 
$. Recall that, for the same reasons, the slope of the conductance-voltage characteristic of 
junction DESVIB is positive at these bias voltages. In addition, at the opening of the second 
inelastic channel, e<£> ~ 2Qi, the reduction of conductance due to vibrationally induced de- 
coherence may also prevail over the increase of conductance due to the onset of non-resonant 
excitation processes. This leads to a drop in the conductance of junction CONVIB at these 
bias voltages. Such a decrease of conduct ance (due to the opening of inelastic transport 



channels) has been reported before 47J, l98HlO(| but only for highly conducting molecular 



junctions (d//d$ > 20/iA/V). We thus conclude that quantum interference effects and vi- 
brationally induced decoherence have a strong impact on the inelastic tunneling spectra of 
a molecular junction. 

Note that antiresonances such as they appear in model CON have been studied extensively 



before and have already been considered in the context of molecular spintronics 87 . 



thermoelectric devices 



or 



85 



86j . Thereby, it is desirable that the antiresonance is located near 



the Fermi energy of the junction. Our results suggest that vibrationally induced decoherence 
puts a strong constraint on these applications. As dynamical excitation and deexcitation 
processes are weakening such antiresonances in both the resonant and the non-resonant 
transport regime, they should be located in a rather narrow range of energies around the 
Fermi energy of the junction, [ep — fii ow ,eF + ^iow], which is determined by the frequency 
Q\ ow of the low-frequency modes of the respective system. 

This concludes our study of basic decoherence mechanisms in vibrationally coupled elec- 
tron transport through single-molecule junctions. For the sake of simplicity, we have so 
far considered a very specific electronic- vibrational coupling scenario (An = 0). In real 
molecular junctions, the intra-state electronic-vibrational coupling will not vanish, except 
if special symmetries are present. The existence of decoherence effects, however, only re- 
quires An 7^ A12, a condition that is fulfilled in many molecular junctions, because typically 
electrons couple state-specific to vibrational modes. 
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D. Suppression of Local Heating by Destructive Quantum Interference Effects 



In Sees. IIII Al - IIII Cl we have shown that destructive interference effects suppress the elec- 
trical current flowing through a molecular conductor. In this section, we demonstrate that 
the respective level of vibrational excitation may also be strongly suppressed, in particular 
in the resonant transport regime, where the level of current-induced v ibrationa l exci tation 



51 
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93 
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94 



126 



127, 



1531 ). 



can be very high (cf. Figs. [7b and [TUb or, e.g., Refs. 

To understand the vibrational excitation characteristics of a molecular conductor, one 
has to take into account excitation and deexcitation processes due to electron transport 
(cf. Figs. - EJi for resonant transport processes and Figs. [T3"r - [T3"b for non-resonant 
transport processes) but also deexcita tion processes, which are associated with electron-hole 



pair creation processes in the leads [154J. Examples for such electron- hole pair creation 
processes are depicted in Fig. [T3b and Fig. [151 which show non-resonant and resonant pair 



creation processes, respectively. Thereby, electrons are not transferred from one electrode 
to the other, as in a transport process, but transfer onto the molecular bridge and, from 
there, back to the electrode where they came from. Although these processes do not directly 
contribute to the electrical current that is flowing through the junction, they re present a 
very important cooling mechanism for the vibrational degrees of freedom 57], [93|, [94J, Il28 |. 
Therefore, as the efficiency of transport processes is, inter alia, determined by the level 
of vibrational excitation, they also have an indirect influence on the electrical transport 

n 

properties of a molecular conductor, which, nonetheless, can be substantial [57J . 

In the present context, it is interesting to note that in model DESVIB electron- hole pair 
creation processes are not suppressed by destructive quantum interference effects, which 
originate from the different signs in the molecule-lead coupling of the two electronic states to 
the right electrode. This is because inelastic processes are, per se, not influenced by quantum 
interference effects due to the specific electronic- vibrational coupling scenario in this junction 
(i.e. because An = 0). As a result, the respective level of vibrational excitation is very 
similar to that of a single vibrationally coupled electronic state. We, therefore, consider 
another model system in this section, model DESVIB2. It is very similar to model DESVIB 
but, in contrast, includes electronic-vibrational coupling also in state 1: An ~ A21 but 
An 7^ A21 (see Tab. [J for a complete list of parameters). Except for molecular juntions with 
special symmteries, this coupling scheme represents a more realistic scenario, in particular 
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FIG. 15: Examples of resonant electron-hole pair creation processes. In contrast to transport 
processes, which involve both electrodes, these processes involve just one of the leads. They are thus 
not subject to the destructive quantum interference effects that are active in model DESVIB and, 
consequently, if destructive interference effects are very pronounced, lead to a strong suppression of 
the current-induced level of vibrational excitation in such a system (e.g. via the resonant excitation 
processes depicted in Figs. Eh and Eh). 

for quasi-degenrate electronic states. As a consequence of the coupling in both states, the 
interaction of the transferred electrons with the vibrational degree of freedom provides less 
accurate which-path information, as inelastic processes can occur in model DESVIB2 also 
with respect to state 1 but with a slightly smaller probability. Inelastic transport processes 
are, therefore, suppressed in junction DESVIB2 by destructive quantum interference effects, 
including excitation and deexcitation processes. In contrast, electron- hole pair creation 
processes are not suppressed, because they involve the coupling of an electronic state to 
just one of the electrodes and not to both, and, even if the respective coupling strength is 
negative, its absolute value is not. 

This has a strong impact on the excitation levels of junction DESVIB2. The correspond- 
ing vibrational excitation characteristic is shown in Fig. [16]by the solid black line. It exhibits 
very similar features as the one shown by model DESVIB but approximately only half of 
the corresponding vibrational excitation levels (cf. Fig. [7b). This is because electron- hole 
pair creation processes represent the dominant class of processes in this system. At low 
temperatures (/cb^l+r << ^i)> they involve only deexcitation processes (as the respective 
excitation processes are Pauli-blocked) while transport processes involve both heating and 
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FIG. 16: (Color online) Vibrational excitation characteristics of the linear molecular conductor 
DESVIB for different values of the polaron-shifted energy level of state 2, ?2- 

cooling processes. Thus, the ratio between heating and cooling processes, w hich is typically 



128 



is shifted 



between 1:2 and 1:1 for resonant electron transport through a single state 
towards much lower values, resulting in much lower levels of vibrational excitation. 

This cooling effect is even more pronounced if the level spacing between the two electronic 
states is smaller. The solid turquoise and red line in Fig. [161 f° r example, show the excitation 
characteristics of model DESVIB2 for smaller splittings of the (polaron-shifted) energies. 
The corresponding excitation levels are substantially smaller. The only major contribution 
that remains in the limit e± — > I2 is the one due to the formation of a polaronic state, 
2~^ mn (^miA„i/fii) (clnCmclcn)-^ (cf. Eq. f fT9|) ). This, however, is only a static contribution 
that originates from charging processes and cannot be used in deexcitation processes by 
other electrons. In other words, the current-induced excitation of the vibrational mode is 
almost completely suppressed in this limit. 

Having in a mind that state 1 and state 2 of model DESVIB2 may represent just sym- 
metric and antisymmetric combinations of localized molecular orbitals (see Figs. and 
[Tb), their electronic structure may be very similar and, accordingly, also the respective 
electronic- vibrational coupling strengths. Indeed, such p airs of electr onic states are often 
found in realistic models of single-molecule junctions [l4, 52, 121 . 122 1. The cooling mech- 



anism described in this section can, therefore, be expected to occur in a large variety of 
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single-molecule junctions. 



E. Decoherence-Induced Temperature Dependence of the Current 

The considerations in Sec. IIII Dl have shown that electron-hole pair creation processes 



facilitate an important coo 



molecule contact 



57|, 



93 



94 



i ng m echanism for the vibrational degrees of freedom of a single- 



128]. In the presence of strong destructive quantum interference 
effects, they may represent even the dominant cooling mechanism. Within the local repre- 
sentation of the linear molecular conductor by two localized molecular orbitals (as in Fig. 
[Tp), this can be understood in an alternative way: If the intramolecular coupling A becomes 
very weak, A — > 0, the left and the right part of the molecular conductor decouple. In that 
sense, the two parts of the molecule behave as being part of the left or right electrode, 
repsectively. Consequently, they acquire the same equilibrium state, that is they thermalize 
with the respective electrode. This thermalization is facilitated by electron-hole pair creation 
processes and is comparable to the thermalization of a molecule adsorbed on a surface with 
the substrate. As a result, the temperature in the electrodes can be used in that regime 
to control the excitation levels of the vibrational modes of the left and the right part of 
the molecular conductor. Consequently, as higher levels of vibrational excitation enhance 
the effect of vibrationally induced decoherence, the temperature in the electrodes can also 
be used to control quantum interference effects in single-molecule junctions. This will be 
demonstrated in this section. 

In general, a variation of the temperature of the electrodes is accompanied by a change 
in the contact geometry of a single-molecule contact, either due to thermal expansion of the 
electrodes or irreversible drifts of single gold atoms. Often, the stability of the contact ge- 
ometry can only be maintained for low temperature variations [ljj]. In order to demonstrate 
control of quantum interference effects via the temperature of the electrodes, we therefore 
consider a model of a single-molecule contact that includes low-frequency vibrational modes, 
which, in the given range of temperatures, exhibit significant changes in their excitation lev- 
els. To this end, we consider a mode with a frequency of Q± = 5meV, which represents a 
frequenc y va l ue a t the lower end of the vibrational spectrum of a larger molecule in a junc- 



tion 



52 



121 



122| . Moreover, we use a range from 10 K to 100 K for the temperature of the 



electrodes in the following, where the corresponding excitation level of the vibrational mode 
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varies between iV v ib = 0.003 and iV v ib = 1.3, respectively. This model system is referred to as 
model DESCOOL (see Tab. [J for a complete list of parameters). Besides the low- frequency 
mode, it includes two electronic states that are almost degenerate, that is I2 — ~£\ <C T. 

The current-temperature characteristic of this model system and the corresponding level 
of vibrational excitation is shown in Fig. [TTJ Thereby, we have used a bias voltage e$ = 
e0.2 V^> 2ej/2, which corresponds to the resonant transport regime of this junction. Note 
that, at this bias voltage, the current level of junction DESCOOL is influence d neither by the 



thermal broadening of the Fermi distribution functions in the leads [145j-[l48| nor vibrational 
sidepeaks that are associated with the onset of resonant excitation processes (Fig. [Sh.). As 
can be seen in Fig. [T7J however, the current and the excitation level of the vibrational 
mode of junction DESCOOL increases substantially with the temperature in the electrodes. 
We attribute this behavior to the quenching of destructive quantum interference effects by 
inelastic transport processes, which is enhanced, via the thermalization mechanism described 
above, for higher temperatures of the electrodes. This shows that the temperature in the 
electrodes can be used to control destructive quantum interference effects in a single-molecule 
junction. Very recently, this control mechanism has been demonstrated by Ballmann et al. 
in mechanically controlled break junction experiments for a variety of different molecules 

Q. 

It is noted that the current levels of model DESCOOl are rather low (tens of pA). This 
is due to the specific choice of the model parameters. To ensure the validity of the effec- 
tive factorization of the single-particle Green's function, Eq. fl2]), we have chosen the level 
broadening V = 0.8 meV significantly smaller than the frequency of the vibrational mode 
Qi = 5meV. While this is a prerequisite for our nonequilibrium Green's function formal- 
ism, the underlying physical processes described in this section are effective also for larger 
molecule-lead coupling strengths. 

IV. CONCLUSION 

We have investigated the role of vibrations in electron transport through single-molecule 
junctions, which are governed by strong quantum interference effects. To this end, we 
have employed both analytical and numerical results that are obtained from nonequilibrium 
Green's function theory 51, 91H9 
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FIG. 17: (Color online) Panel (a): Current-temperature characteristic of the linear molecular con- 
ductor DESCOOL. Panel (b): Vibrational excitation characteristic corresponding to the current- 
temperature characteristic shown in Panel (a). 

Our results show that electronic-vibrational coupling has a strong and non-trivial in- 
fluence on the transport properties of a molecular conductor, in particular, if quantum 
interference effects play a dominant role. We have investigated molecular junctions, where 
pronounced quantum interference effects arise due to the existence of quasidegenerate states, 
which provide different pathways for electron transport. Our results show that electronic- 



40 



vibrational coupling may strongly influence quantum interference effects in molecular junc- 
tions. This is because the electronic states of a molecule couple, in general, very specifically 
to the respective vibrational degrees of freedom. The different coupling provides information 
about the pathway taken by the electron. This which-path information quenches quantum 
interference effects. This decoherence mechanism is more pronounced for higher levels of vi- 
brational excitation or, equivalently, higher effective temperatures of the vibrations, because 
the probability for inelastic transport processes increases with temperature. 

At high bias voltages, where resonant transport processes may result in high levels of 
vibrational excitation, vibrationally induced decoherence is therefore very pronounced and 
quantum interference effects play only a minor role. At the onset of the resonant transport 
regime, they may result in sizeable effects, such as, for example, a reorganization of step 
heights in the current- voltage characteristics or a modification of the widths of these steps. 
At low bias voltages in the non-resonant transport regime, where the excitation levels of 
the vibrational modes are typically rather low, vibrationally induced decoherence is less 
pronounced. In this regime, interference effects may, nevertheless, result in a strong modi- 
fication of the signals that are associated with inelastic electron tunneling. For example, in 
the presence of destructive interference effects, the onset of inelastic processes is, in general, 
more pronounced while in the presence of constructive interference effects the respective sig- 
nals are less pronounced or even reversed, that is, negative jumps in the conductance- voltage 
characteristic appear, even if th e conductance of the molecule is much smaller than the con- 



ductance quantum 47|, l98Ml00| . A correct and thorough analysis of vibrational signals in 



the transport characteristics of a single-molecule contact needs to take into account these 
effects, which may lead to strong deviations from the commonly employed Franck-Condon 
picture. 

We have also elucidated the role of electron-hole pair creation processes in the pres- 
ence of strong destructive quantum interference effects. In contrast to transport processes, 
electron-hole pair creation processes are not suppressed by interference effects and, therefore, 
constitute the dominant inelastic processes. At low temperatures of the electrodes, they are 
solely cooling the vibrational degrees of freedom. In this case, the effective temperature of 
the vibrational modes is only very weakly influenced by inelastic electron transport processes 



and is rather determined by the temperature in the electrodes 



60( ■ Because the level of vi- 



brational excitation determines, inter alia, the efficiency of transport processes and, thus, 
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also the effect of vibrationally induced decoherence in such a junction, electron-hole pair 
creation processes thus facilitate an effective mechanism to control quantum interference 
effects in single-molecule junctions by the temperature in the electrodes. This has recently 
been demonstrated in experiments on a variety of different molecular junctions by Ballmann 
et al. ljj. 

While we have presented a comprehensive study of vibrationally induced decoherence 
in single- molecule junctions, further research is required that addresses not only the anti- 
adiabatic regime but also the adiabatic and the respective cross-over regime. Moreover, 
because quantum interference effects appear in particular for quasi-degenerate electronic 



states, it is a 
coupling 



18 



so intere sting to investigate the effect of non-adiabatic electronic-vibrational 
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143| in these systems. Thereby, theoretical methods need to be applied 



or developed that allow a nonperturbative description of electronic-vibrational coupling as 
well as higher order effects. 
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Appendix A: Transmission function for junction DESVIB 

In this appendix, we derive the transmission function of model DESVIB and the corre- 
sponding incoherent and interference term (cf. Eqs. ( 1381) . Sec. IIIIBl) . To this end, we employ 
the wide-band approximation, T « r L) ii(// L ), and consider large bias voltages ($ — > 00). 
Furthermore, we evaluate the vibrational degree of freedom in thermal equilibrium, that is 
we restrict it to an effective (finite) temperature T. Given these assumptions, the real-time 
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projections of the corresponding electronic self-energy matrices (11 II) can be written as 

, iT iTA , 

SfJ e = , (Al) 

iAT iT 

£>(e) = 0, (A2) 
£<(e) = 0, (A3) 
, -iT iTA \ 

XZ(e) = \ , (A4) 

iAT -iT J 

, —iT — iTA i 

S r L 6 = =(S£(6)) f , (A5) 

— iT iTA 

S r R W = I I =(SM) t , (A6) 



zTtl -zr 



with 



A = e -^T (2iW1) , (A7) 

iVvib = t • (A8) 

e fe B T — 1 

The retarded/advanced projection of the electronic part of the single-particle Green's func- 
tion (cf. Eq. ffl2|) ) is therefore given by 

£r/a (e) = ( (<0 \ / \ 

V atie)) \ o 

and, according to the Keldysh equations (fl3l) . we write the corresponding greater Green's 
function as 

-iT\g\{e)\ 2 iTAg\{e)gl{e) 
iTAgl{e)gt{e) -iT\gl{e)\ 2 

The greater projection of the single-particle Green's function matrix (see Eq. (Q) thus reads 

-iT\gUe)\ 2 iTA 2 gUe)g%(e) 
G > (e) = mK M yiK )y2K ' , (A10 

l TA 2 gl{e)gt{e) -iT A 2 ^T=-oo Bi\ 9 r 2 (e - IQ 1 )\ 2 
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with 



B l = h (^2^^N vih (N vlh + l)^ e^T (All) 
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and /; the lih Bessel function of the first 
of the shift operator correlation function 

^A 2 



rind . Thereby, we have used the lesser projection 
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(X m (t)X* n {lf))> = 8 m2 5, : 



l 2exp 



^ (i\U(e^') -!) + (! + N vih )(e-^W - 1)) 



(A12) 
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and convoluted this function with the lesser projection of the electronic part of the single- 
particle Green's function according to Eq. fl9]). Finally, the transmission function is obtained 
using the Meir-Wingreen-like formula ( 12 ip . which can be rewritten as 



(A13) 



with 
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(A14) 
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For large bias voltages, this expression can be used to define a transmission function by 



de 
2^ 



Tr 



Ef'<(6)G>(e) 



= 2e 



de 
2^ 



*(e). 



This transmission function is given by 
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= Tr 
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The respective incoherent term is obtained by the neglecting the off-diagonal elements of 
El m Eq. (1A17I) (cf. the discussion given in Sec. IHI Aj) 



t iac (e) = r 2 |^( e )| 2 + r 2 A 2 Bi\gl(e-in 



IJI 2 ) 



l=— oo 



which also allows to identify the corresponding interference term as 

W e ) = t(e) - tinc(e) 

= -T 2 A 2 (gl(e)gl(e)+gl(e)g*(e)). 



(A18) 



(A19) 
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Appendix B: Transmission function for junction CONVIB 



In this appendix, we derive the transmission function of model CONVIB and the corre- 
sponding incoherent and interference term (cf. Eqs. Sec. IIIIBl) . To this end, we employ 
the same approximations as in appendix The corresponding real-time projections of the 
electronic self-energy matrices (TTTT) are given by 
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leading to the following retarded/advanced projection of the electronic part of the single- 
particle Green's function (cf. Eq. f|T2|) ): 
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According to the Keldysh equations ( fl3l) . we can write the greater projection of the (elec- 
tronic) Green's function as 
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-e l + tr){e-e 2 + tT) + r 2 A 2 


Zn{e) 




e — 


e 2 + z r| 2 -r 2 A 2 , 


Z 12 (e) 


= A(e 


-e 1 ){E-e 2 )-T 2 A{l-A 2 ), 


Z 21 (e) 


= A(e 


-^(e-^-P^l-A 2 ), 


Z22 (e) 




e — 


e x + iT\ 2 -T 2 A 2 . 



(B8) 
(B9) 
(BIO) 
(Bll) 
(B12) 

Convolution of this expression with the greater projection of the shift operator correlation 
functions flA12h leads, according to Eq. 0, to the greater projection of the single-particle 
Greensfunction matrix 



G>(e) = -iT 



c(e)Z n (e) 



Ac{e)Z 12 {e) 



Ac(e)Z 21 (e) A 2 J2Z-oo B ^ ~ Mx)Z 22 {e - IQ X ) 
The corresponding transmission function, as defined by Eq. (IA16[) . thus reads 



(B13) 



t(e) = Tr 



< Q )><(,\r<> 



T 2 Tr 



e)G>(e) j 

1 1 W c(e)Z n (e) Ac(e)Z 12 (e) 
1 1 I \ Ac(e)Z 21 (e) A 2 ^ c ( e ~~ l^i)Z 22 (e - lQ t ) 



(B14) 



T 2 c{e)Z u {e) + Ac(e)Z 21 (e) + Ac(e)Z 12 (e) + A 2 ^ c ( e " ^i)^ 2 (e - ^1) 



i=— 00 



Considering the contributions of the off-diagonal elements of £l' < -°' ) in Eq. (1B14I) (cf. the 
discussion given in Sec. IHI Al) . it can be decomposed into the interference term 



£ int (e) = AT 2 c(e)(Z 21 (e) + Z 12 (e)), 

and the incoherent term 

/ 00 
t inc (e) = r 2 lc(e)Z ll (e) + A 2 ^ B lC (e - IQ^Z^e - IQ, 



(B15) 



(B16) 



i=— 00 
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